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The  Fourth  periodical  report  concentrated  on  the 
numerical  solution  of  H.L.S.  equation  with  varying  coeff lclenta. 

The  nuaerlcal  treatment  la  baaed  on  the  analytical 
solution  developed  In  the  third  periodical  report.  The  solution 
la  valid  using  periodic  boundaries  and  a  very  alow  variation  of 
the  bottoa  compared  with  the  period's  length.  The  solution  was 
compared  with  other  nuaerlcal  schemes  and  yield  satisfactory 
results. 

All  technical  details  are  presented  In  the 
following  paper. 
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zvoumoM  ot  an  un stable  hater  have  train  moving  over,  a  slowly 

VARYING  BOTTOM 
(Part  II) 

by  t.  I us  la  and  Y.  Frost lg. 


Introduction 

- - Tha  ala  of  this  pcpar  Is  to  present  quantitative  Information 

about  the  approximate  solution  of  tha  non-linear  Schroedlnger 
-  (M. L.-frrfr^equatlon  with  varying  coefficients  and  periodic  boundary 


conditions 


►.^coeparlng  It 


srical  solutions.' 


In  section  I  we  present  |  brief  review  of  the  aethod^ developed 
is  <2>,  whirls  based  la  tbs  evaluation  of  an  ^Initial  disturbance 
at  every  point,  and  local  application  of  the  analytical  solution . 
glveu  by  K.  Stlassele  at  el.  <t>.  In  section  It -wo  provldd*  quantitative 
intonation  about  the  variation  of  the  local  ^initial  disturbance*^* 
and  in  a  or  tine  fl^we  present  the  comparison  of  the  solution 
dswelopad  hHft  with  nunarical  solutions. _ _ _ 


In  the  Previous 


l.  Key lew  of  the  Math 


The  problan  to  be  solved  la: 

‘♦x  +  *n  +  “  0 

with  tha  boundary  condition 

*1  -  1  +■  2*eUcoa(2*T) 

l*-0 


2- 


Pol loving  Stlaaanla  and  bMijmkjv  (1)  we  approximate  the  eolutlon 
of  (1.1),  (1.2)  by  the  truncated  Fourier  aeriea: 

t*(T,X)  -ID  (X)e2,lnT  (1.3) 

n-l  “ 

where  D  (X)  la  the  complex  Fourier  coefficient  of  aodulua  Id  | 
n  d 

and  argunant  an: 


The  boundary  condition  at  X  -  0  lapliaa  that 
|Bo(o)|  -  1,  aQ(o)  -  0 

IDjCo) I  -  5,  el(o)  -  a 

Fran  eubetltutloa  ef  (1.3)  Into  the  N.L.S.  (1.1)  the 
following  eat  of  equatlooa  la  obtained 


(1.4) 


(1.5) 


dD 

1  +  M(X)  [(|dJ2  +  *|Dj|2)Do  ♦  2DjD*]  -  0 

(1.6) 

dD±. 

i  -5T+  w(X)  [«|D  |2  +  3|D  I2  -  WDj  ♦  D2D*]  -  0 


where  P  -  8«2/u  .  The  range  of  variation  of  P  la  1  <  P  <  6.  (aee 
Stlaaanla,  at  al.  (1)). 

Aeauaing  that  the  bo t toe  varlea  eo  (lowly  that  It  nay 
be  conaldared  "locally  conatant”,  then  the  eolutlon  given  by 
Stlaaanla  at  al.  (1)  for  conatant  depth  can  be  applied  locally. 


(1.16) 


c  -  j  [4(l+282)-P]  {  1-  [1-  3%> 

Jj 

d  -  p  {  i  -  [  i  +  :>  a.i7) 

p 

If  I(P)  >  0  or  vice-versa,  d  given  by  (1.16)  end  c  by  (1.17)  If 
I(P)  <  0. 

The  approxlnats  solution  depends  on  two  psraneters:  I(P)  and  the 
local  initial  condition  for  Z:  E(P).  Given  1(P)  sad  *(P) ,  Z(X)  Is 
determined  from  (1.8)  to  (1.17).  At  the  case  of  constant  depth: 


£(P)  a  const  •  282 


(1.18) 


and  I(P)  “  const  •  62[  2$2  +  4(1  +  cos  2a)  -  2P] 


(1.19) 


Voder  the  assumption  that  the  bottom  Is  "locally  constant" 

In  the  sense  previously  mentioned,  hs  derived  In  (2)  the  following 
approximate  differential  equation  for  I(P). 


1^ 

2  ln( - T  -  -  --  ) 


»p  y7 


•  ^ 


«  2PZ(»-P)2 

1b  (2P-0l<P) 


(1.20) 


with  the  Initial  condition  (1.19). 


I 
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Equation  (1.20)  can  ba  aaally  solved  numerically. 

In  order  to  derive  equation  (1.20)  we  assuned  in  that  the  initial 
condition  of  Z(X)  for  local  disturbance  satisfies  the  condition 

e  -  e(P)  «  1.  We  have  chosen  as  the  local  Initial  condition 
e(P),  the  minimal  positive  value  of  the  oscllatory  function  Z(X), 

The  minimal  value  of  Z(X)  la  found  by  equalling  to  aero  Its 
derivative  given  by 

H  "  -^2ZU+2B?-Z)]2-[  |z2+[P-2(1+2B2)]Z  -  I(P)]2  (1.21) 

(see  Stlassnie  et  al.  (1)). 

If  I(P)  >  0  the  minimal  positive  value  of  Z  is  given  by 

e(P)  -  }  [*(l+2B2)-P]  {!-  [  l - kW)  ]V=c  (1.22) 

*(1+2B2)-P 

and  from  (1.9)  It  can  be  seen  that  it  corresponds  to  cos2(oo-ol>  «  1 
If  I(P)  <  0  the  minimal  positive  value  of  Z  is  given  by 

e(P>  -  P  {  1  -  [  1  +  }  =  c  (1.23) 

P 

and  it  corresponds  to  cos  2  (a  -a.)  •  -1. 

O  1 

The  approximate  wave  envelope*  j(T,X)  is  given  by 
2  2 

|*  |  -  (-2I(P)  +  2[*(1+2BZ)-P]Z  -  7Z2  + 

+  [S  (2I(P)  +  2PZ-Z2)*1  +  *Zcos2»T]2}  /8Z. 

where  8  is  the  sign  of  cos  (a.-o  ). 

I  o 


(1.2*) 


XI.  The  behavior  of  I(P) 

The  right  hand  aide  of  equation  (1.20)  la  always  negative  for  the 
range  of  values  of  P  conslderated  1  <  P  <  A,  yielding  that  I(P)  la  a 
decreasing  function  of  P.  On  the  other  hand,  It  can  be  easily  seen  that 
I(P)  SO  is  a  solution  of  equation  (1.20).  That  means  that  if  the  Initial 
condition  X(PQ)  is  positive,  the  solution  remains  positive,  and 
aslmptotlc  to  the  horizontal  axis.  If  the  Initial  condition  Is 
negative,  the  solution  remains  negative.  The  expected  behavior  of  the  solution 
is  given  in  figure.  1. 


Figure  1 
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I f  I(P)  -  0  then  eqs.  (1.16)  and  (1.17)  yield  c  =  d  =  0  and 
from  (J.13)  we  obtain  k  5  1.  Then  the  period  of  the  elliptic 
functions  in  (1.11)  la  Infinite  (see  (3)). 

The  solution  Z (X)  in  this  case  is  given  by  (1.10)  where 


i  -  (®±.  .£)  «gh  v 
^  va-c  b^  sschzy 

i  +  (i±.  £  *&h2y_ 
'a-c  b  .2 
sech  y 


where  a  and  b  are  given  by  (1.14)  and  (1.15)  respectively.  From  (1.10)  and 
(2.1)  it  can  be  seen  that  Z(X)-0(e).  This  solution  is  unstable  to 
disturbances  in  the  parameter  I(P),  giving  rise  to  the  oscillatory 
solution  (1.10,(1.11)  which  grows  from  its  minimal  value  Z  -  c  •  0(E), 


to  its  maximal  value  z  -  b. 


For  this  reason  all  the  numerical  calculations  are  unstable 
near  the  value  I(P)  «  0.  On  the  other  band,  when  I(P)  tends  to  zero, 
the  period  of  the  oscillatory  solution  Z(X)  tends  to  infinity,  and 
then  our  assumption  that  the  bottom  is  nearly  constant  in  a  period 
of  the  solution  does  not  hold  any  more. 

These  reasons  Impose  new  restrictions  in  F.  The  permissible  range 
of  variation  of  P  must  ensure  the  conditions  |l(P)|  »  1,  but  1(P) 
not  nearly  zero. 

From  the  numerical  results  we  actually  Impose  the  condition 

|i(p)J  >o(io"3). 

The  solution  of  equation  (1.20)  is  compared  with  the  one 


obtained  by  solving  numerically  the  set  (1.6)  and  using  the 
equality 


I(P)  -  |D1j2-[2|Dl|2+4|B0i2-2P44(Do|2co.2(a1-ao)] 


(2.2) 


Equation  (1.16)  traa  aolvad  by  neana  of  the  second  order  Range 
Kutta  net hod,  and  by  the  trapezoidal  rule,  obtaining  identical  results. 
System  (1.6)  was  solved  by  the  trapezoidal  rule.  In  both  cases 
we  considered  the  simple  case  where  P  is  linear  on  X: 

P  «  ?o  +  0.1X;  Pq  -  2  starting  at  X  ■  0.  The  results  are  presented 
In  figure  2.  Full  lines  represent  the  solutions  obtained 
by  (1.6)  and  (2.2),  while  dotted  lines  are  the  results  from 
numerical  solutions  of  equation  (1.20). 
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III.  Comparison  between  the  solution  presented  In  section  I 
and  numerical  solutions 

In  order  to  appraise  the  validity  of  the  behavior  predicted 

by  equation  (1.20)  and  the  analytical  solution  (1.7)  to  (1.17),  we 

compare  It  with  a  reference  solution  obtained  by  means  of  numerical 

schemes.  In  the  first  place  we  compare  it  with  the  numerical 

solution  of  the  set  (1.6)  obtained  by  means  of  the  trapezoidal 

rule.  Figures  3  and  4  show  the  wave  envelope  |i|>(0,X)|  for  the  case 

P  -  P  +  0.1X,  p  -  2,  a  -  0,  S  “  0.1. 
o  o 

These  values  correspond  to  the  Initial  condition  I(PQ)  *  0.0402. 

Full  lines  represent  the  solution  of  the  set  (1.6)  while  dotted  lines 

are  the  results  obtained  by  solving  equation  (1.20)  and  applying 

(1.22),  (1.11)  to  (1.17)  and  (1.24)  at  P  -  2,  P  -  0.5  and  P  -  2.16. 

Figures  5  and  6  represent  the  wave  envelope  |i|((0,T)j  when  a  •  nil, 

B  -  0.1  P  -  P  40. IX,  P  -  2,  (1(P  )  -  -0.0998).  The  full  lines 
0  0  o 

represent  the  solution  of  set  (1.6),  and  the  dotted  ones  the 
solution  given  by  (1.20),  (1.23),  (1.10)  to  (1.17)  and  (1.24)  at 
P  -  2  and  P  -  2.5.  In  both  cases  a  very  satisfactory  correspondence 
of  the  wave  envelop 's  shape  is  observed. 

On  the  other  hand,  the  solution  of  set  (1.6)  is  compared 
to  the  numerical  solution  of  the  N.L.S.  equation  when  P  •  2  +  0.1X. 

In  order  to  aolve  the  N.L.S.  equation  two  Independent  alternative 
numerical  schemes  were  employed,  obtaining  Identical  results. 

The  first  one  is  the  Crank  Nicholson  scheme 


_ki_ 
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♦r1  -  *“  *t\  -  <i + cl  <-i  - 2*: ♦ Cv 

i 1  +  *i  (A1 - ^ - =*■  +  -*-1 - ^ — e*-) 


AX 


(AT) 


,,n+l  .  ,n. 

(*,-  +  *«) 


(AT) 


(3.1) 


where  *"  -  *(Tj,Xn),  u“  -  utt^),  Tj  »  (j-l)AT,  Xr  -  nAX,  AT  are 

J-l  equal  segments  chat  span  the  Interval  0  <  T  <  !s  and  AX  are  equal 

Intervals  that  span  the  X  coordinate.  The  scheme  Is  subject  to  the 

Initial  condition  (1.2)  and  to  the  boundary  conditions  ijin  -  >|n ; 

o  2 

*n  -  jn  . 
vj+i  j-r 

In  the  second  approach  wo  seek  an  approximate  solution  of  the 

form 


/(T.X) 


N 

I  D  (X)e 
«-(N-l)  n 

Jf, 


2irinT 


(3.2) 


Given  the  set  f  (Tj,X)  j  •  -(N-l) , . . . ,0,1, . . .N,  the  corresponding 

Fourier  coef f lclentsDOx^^  , . . . ,  DU)  ' . (x)  are  found  by  means 

of  the  Fast  Fourier  Transform,  and  conversely,  given  -1)  M . 

Dq  , ..  ,Djj  (x)  the  set  <|N(X  ,Xj=-(N-l) ,. .  ,0,1, . .  ,N  is  found  by  means 
of  the  Inverse  Fourier  Transform. 

Substitution  of  (3.2)  In  the  N.L.S.  equation  (1.1)  yields: 


IdD 

M  +  Kin)4®  -  iuA  ■  0 
®  n  n 


(3.3) 


D  -  D 
-n  n 


(n  -  0,1,. ...N) 


where  A  »  A  (X),  n  «  -(N-l)  ,. .  ,0,. .  ,N  are  the  Fourier  coefficients 

n  n 

of  the  set 

|/(Tj,X)|2.*N(Tj,X),  j  -  -(N-l)  ,. .  ,0, .  .N  . 


i 


The  solution  of  (3.1)  Is  given  by 


D.  (x)  -  [D  (0)  +  ip  /*A  (Oo1(2lra)2tdt3  e"1(2,1)  *  (3.4) 

-n  o 

which  can  be  rewritten  in  the  fons 

D*  (x+AX)  -  [D  (x)  +  ipe‘1(2TO)2x  1***  e1(2TO)2tA  (t)dt]  e'112'^ 
xn  n  n 

n  “  0,1, ... ,N  (3.5) 

Assuming  that  the  chosen  increnent  AX  la  sufficiently  snail,  the 
Integral  in  (3.5)  is  conputed  approximately  as 
x+Ax  .2 

»i[A  (x)  +  A  (x+AX)]  /  e1(Zlrn)  *dt  (3.6) 

“  n  , 

Substitution  of  (3.6)  into  (3.5)  leads  to  the  numerical  scheme 


D  (X+AX)  -  D  (X)  +  AXip(A  (X+AX)  +  A  (Xj/2 
o  o  o  o 

D*.  (X+AX)  -  f  0  (X)  +  (1-f  )(A  (X+AX)+A  (X))/(P-n2) 
*n  n  n  an  n 

n  -  1,2 . « 

,  ,  -i(2xn)2AX 

wh«r«  f  •  e 

n 


(3.7) 


Starting  with  the  known  values  of  D  at  level  X,  we  proceed 

It 

to  c caputs  A  (X)  by  naans  of  the  Fast  Fourier  Transform.  We  use  this 

U 

value  as  a  first  guess  for  A  (X+AX)  and  obtain  an  estimate  for  D  (X+AX) 

n  n 

fren  (3.7).  The  Fast  Fourier  Transform  and  expression  (3.7)  are  then 


applied  Iteratively  until  no  change  between  two  aucceaaive  estimates 

are  detected.  The  process  la  then  advanced  a  further  step  In  X.  The  first 

step  uses  as  starting  values  the  Initial  condition  (1.2). 

Figures  7  and  8  show  the  solution  of  the  N.L.S .  equation  when 

P  -  P  +  0.1X.  P  »  2,  for  the  initial  conditions  o«  0,  B  *  0.1 
o  o 

(I(P  )»0.0402)  and  a-  0,  8*  0.05  (I(P  >  ”  0.01000125)  respectively, 

O  o 

while  in  figures  9  and  10  it  can  be  seen  the  solution  for  the  initial 
conditions  a  »  s/2,  6  -  0.1  (I(Pq)  “  -0.0998),  Full  lines  in  these 
figures  represent  the  solution  of  the  N.L.S.  equation,  while  dotted  lines 
represent  the  solution  of  system  (1.6). 


